High threshold universal quantum computation on the surface code 
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We present a comprehensive and self-contained simplified review of the quantum computing 
scheme of [3, 01 1 which features a 2-D nearest neighbor coupled lattice of qubits, a threshold error 
rate approaching 1%, natural asymmetric and adjustable strength error correction and low over- 
head, arbitrarily long-range logical gates. These features make it by far the best and most practical 
quantum computing scheme devised to date. We restrict the discussion to direct manipulation of 
the surface code using the stabilizer formalism, both of which we also briefly review, to make the 
scheme accessible to a broad audience. 



I. INTRODUCTION 

Classical computers manipulate bits that can be ex- 
clusively or 1. Quantum computers manipulate quan- 
tum bits (qubits) that can be placed in arbitrary su- 
perpositions a|0) + f3\l) and entangled with one another 
(|00) + |ll})/\/5. This additional flexibility provides both 
additional computing power and additional challenges 
when attempting to correct the now quantum errors in 
the computer. An extremely efficient scheme for quan- 
tum error correction and fault-tolerant quantum compu- 
tation is required to correct these errors without mak- 
ing unphysical demands on the underlying hardware and 
without introducing excessive time overhead and thus 
wasting a significant amount of the potential performance 
increase. 

This paper is a simplified review of the quantum com- 
puting scheme of (T|, [3| ■ The scheme requires a 2-D 
square lattice of nearest neighbor coupled qubits with 
initialization, readout, memory and quantum gates all 
operating with error rates less than approximately 1% 
— the least challenging set of physical requirements de- 
vised to date. Furthermore, despite the modest physical 
requirements, logical qubits (qubits of data distributed 
over many physical qubits and protected by error cor- 
rection) can be interacted over arbitrarily large distances 
with time overhead only growing logarithmically in their 
separation. This is remarkable as most nearest neighbor 
quantum computing schemes are associated with a time 
overhead that grows linearly with logical qubit separa- 
tion. Finally, in most, if not all, physical quantum com- 
puter technologies, bit-flips |0) «-> |1) are less likely than 
phase-flips |1) — |1) opening the door for asymmetric 
error correction schemes that make use of fewer physical 
qubits to preserve a given amount of quantum data with 
a given confidence level. Practically, it is also helpful if 
additional error correction resources can be dynamically 
allocated to critical data during the quantum computa- 
tion. The scheme we review permits both asymmetric 
and dynamic error correction in a natural manner. 

A number of technologies are well-suited to implement- 
ing surface code quantum computing. Proposals exist 



for 2-D arrays of qubits making use of superconductors 
[H, |i[ and semiconductor nanophotonics Q . An equiva- 
lent measurement based version of the scheme calling for 
a 3-D cluster state fl L 1, @ could be implemented using 
photonic modules U 8 ] or 10n t ra P s 0] • 

The discussion is organized as follows. In Section |TT] 
we briefly review the stabilizer formalism of quantum 
computing 0. Section [m] briefly reviews the surface 
code fill ], which forms the error correction substrate of 
everything that follows. Logical qubits are introduced 
into the surface code in Section HVi along with their ini- 
tialization, measurement, and basic logical operations. 
Section [V] describes logical CNOT in detail. Section IVII 
completes the universal set of logical gates with a discus- 
sion of state injection, state distillation and appropriate 
quantum circuits making use of the distilled states. An 
efficient implementation of logical Hadamard inspired by 
[l2| that avoids the extensive machinery of Section IVII 
is described in Section fVIII Section [Villi then describes 
simulations used to estimate the threshold error rates 
of physical qubit initialization, measurement, memory 
and two-qubit gates. Looking further into the future, 
Section IIXI discusses distributed quantum computing to 
make it clear that impractically large 2-D square lattices 
of qubits are not required to tackle problems of interest- 
ing size. Section IXl summarizes the discussion and points 
to further reading. 



II. QUANTUM STATES AND STABILIZERS 

A quantum state can be specified in a number of equiv- 
alent ways. One of the most common is to choose a 
basis and express the state as a state vector such as 
(| 00) + \U))/V2. In this review, it will be much more 
convenient to express this state as the unique simultane- 
ous + 1 eigenvector of the commuting operators X g) X 
and Z®Z. Such operators are called stabilizers [l(| . This 
entire review is based on the manipulation of stabilizers. 

Any set of n mutually commuting and independent 
operators over n qubits has a unique simultaneous +1 
cigenstate. We will restrict our attention to stabilizers 
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that are a tensor product of the identity operator / and 
the Pauli matrices X, Y, Z (with Y = XZ real). A 
set of such stabilizers cannot be used to specify an ar- 
bitrary quantum state, though a sufficiently broad range 
of states can be specified for most of our purposes. See 
Section IVI Al for a simple extension to the stabilizer for- 
malism permitting arbitrary states to be specified. 

Consider a set of stabilizers M specifying state \tp). 
Suppose we wish to apply an operator U to state \ip). 
If we consider U\tp) = UMWU\^)), we can see that the 
new set of stabilizers will be UMW . To give an explicit 
example, 

M = Z1Z2, Z2Z3, Z 3 Z 4l X1X2X3X4 (1) 
U = X 2 (2) 
UMU f = -Z 1 Z 2 ,-Z 2 Z 3 ,Z 3 Z i ,X 1 X 2 X 3 X 4 . (3) 

In addition to unitary manipulation, we will frequently 
discuss measurement of a given operator, for example X, 
Z or some more complicated operator involving a larger 
tensor product. A very simple example is a single qubit 
in an unknown state with stabilizer I and the subsequent 
measurement of the Z operator. We will write the sta- 
bilizer of a qubit after such a measurement as ±,Z. Note 
that the probabilities of the two possible measurement 
outcomes, the +1 and —1 cigcnstatcs of Z, are typically 
not recorded in the stabilizer formalism, just their possi- 
bility. Note also that given any operator there is always 
a nonzero probability of obtaining at least one of the two 
cigenstates. 

Care needs be taken when measuring if other nontrivial 
stabilizers are present. There arc three cases to consider. 
If the operator to be measured can be expressed as a 
product of stabilizers, no change is made to the stabilizers 
as we already have an eigenstate of the operator. For 
example, if we have two qubits and stabilizers Z\ and 
— Z 2 , measuring the Z\Z 2 operator will always give the 
— 1 eigenstate. 

If the operator to be measured cannot be expressed as a 
product of stabilizers and commutes with each stabilizer, 
the operator is added to the list of stabilizers with a sign 
that depends on whether we have projected into the +1 
or —1 eigenstate of the operator. For example, if we have 
three qubits and stabilizers Z\Z 2 and Z 2 Z 3 , measuring 
Z 2 will yield one of the ±1 cigenstates, meaning we will 
introduce the new stabilizer ±Z 2 . Note again that the 
probability of the two outcomes is neither recorded nor 
known. 

Finally, if the operator cannot be expressed as a prod- 
uct of stabilizers and anticommutes with one or more sta- 
bilizers, the second and subsequent anticommuting sta- 
bilizers are multiplied by the first anticommuting stabi- 
lizer to form commuting stabilizers, and the first anti- 
commuting stabilizer is replaced with the operator being 
measured, again with sign depending on which state we 
have projected into. For example, if we have four qubits 
and stabilizers Z\Z 2 , Z 2 Z 3 and Z3Z4, to measure A3 we 
first multiply Z 3 Z 4 by Z 2 Z 3 and then replace Z 2 Z 3 with 
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FIG. 1: Basic layout of surface code data qubits, each repre- 
sented by a circle. A data qubit is located at the center of 
each edge of a square lattice. The square lattice is a guide for 
the eye only, it does not represent interactions. 



±A3 to give the new set of stabilizers ZiZ 2 , ±X 3 and 
Z 2 Z 4 . In this instance we know that the probability of 
the two outcomes is equal as, given any state \ip) stabi- 
lized by some operator S and any operator to measure 
M such that MS = -SM, we have 
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meaning we have an equal superposition of the ±1 eigen- 
statcs of M. 

Many other examples of measurements falling into each 
of these three categories will be discussed in subsequent 
sections. 



III. THE SURFACE CODE 

The surface code was first presented in [ll|. A small 
surface showing the basic layout of qubits, a square grid 
with qubits on each edge, is shown in Fig. Q] The stabi- 
lizers of this surface are 

A0A2, X0X1X3, X1X4, A2X5A7, A3A5A6A8, 

X 4 XqXq, XjXiq, X&X\qX\i, A9A11 

ZqZ 2 Z 3 Zs, Z\Z 3 Z 4 Z§, Z^Z-jZsZiq, Z§Z%Z§Z\\ (5) 

These correspond to a tensor product of Z around each 
face and X around each vertex. Note that XgXu can be 
expressed as a product of the other X stabilizers. This 
leaves 12 independent stabilizers on 12 qubits implying a 
unique state. Given a 11; by ft, face surface, in general there 
will be 2wh + w + h qubits and independent stabilizers. 

Not shown in Fig. [1] are additional qubits on each face 
and vertex that enable one to check the sign of the associ- 
ated stabilizer. These additional syndrome qubits make 
the lattice a simple nearest neighbor connected square 
lattice. Discussion of the syndrome qubits and the quan- 
tum circuits used to extract the signs of the stabilizers 
will be deferred until Section [Villi 

If no errors of any kind occur, the surface remains in 
the simultaneous +1 eigenstate of every stabilizer. When 
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FIG. 2: (Color online). Effect of single bit-flips (X) and 
phase-flips (Z) on the surface code — adjacent face and vertex 
stabilizers are made negative, respectively. 



discussing errors, we will restrict our attention to bit-flips 
and phase- flips. Very general noise can be tolerated with 
just the ability to correct these two types of errors [HI]. 
Fig. [2] shows the effect of single bit-flips and phase-flips 
on the surface code — the adjacent stabilizers become 
negative. If we could reliably detect when a stabilizer 
becomes negative, this clearly would be sufficient to pin- 
point and then correct these single errors. 

Two additional complications need to be accounted for. 
Firstly, it is possible for long chains of errors to occur. 
Secondly, it is possible for the reported eigenvalue of a 
given stabilizer to be wrong. Both of these situations are 
illustrated in Fig. [5] To cope with these complications, 
we keep track of every time the reported eigenvalue of 
each stabilizer changes. Without loss of generality, let 
us focus solely on Z stabilizers, which detect bit-flips, as 
both types of errors are treated independently. 

Fig. [3^ gives an example of appropriate Z stabilizer 
information. In practice, correction is delayed for as long 
as possible and pairs of flipped syndromes are then con- 
nected by paths in space and time or "matched" such 
that the total number of edges used is minimal, as shown 
in Fig. [5}x Polynomial time minimum weight match- 
ing algorithms exist (3], hence this can be done effi- 
ciently. Note that X errors can be matched to smooth 
boundaries and Z errors to rough boundaries of the sur- 
face. A smooth boundary is a boundary with four term 
Z stabilizers and three term X stabilizers as shown in 
Fig. [3J A rough boundary is a boundary with four term 
X stabilizers and three term Z stabilizers, also shown in 
Fig. [3l Given a minimum weight matching, bit-flips are 
applied to the spacclike edges to correct the errors with 
high probability. Further discussion of the details of error 
correction will be delayed until Section IVIIII 

Initialization of the surface code substrate is not com- 
pletely trivial. If every qubit is prepared in the |0) state, 
we automatically have the +1 cigenstate of every Z sta- 
bilizer, but when we measure the X stabilizers the eigen- 
states will be randomly positive and negative. For sim- 



FIG. 3: (Color online). Examples of smooth and rough 
boundaries including a chain of X errors ending in a smooth 
boundary without changing the sign of any stabilizers, and 
a chain of Z errors ending in a rough boundary, also with- 
out leaving any evidence of its presence. X/Z stabilizers are 
represented by green/blue shapes, respectively. 
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FIG. 4: (Color online). Surface code suffering from multi- 
ple errors (indicated by labeled red dots) and an incorrect 
syndrome measurement (indicated by the dashed red circle). 



plicity, we choose to treat the random negative eigenval- 
ues as errors and correct them. 



IV. LOGICAL QUBITS 

Armed with the surface code described in the previous 
section, we can now discuss logical qubits. The simplest 
logical qubit consists of a single face where we stop mea- 
suring the associated Z stabilizer. This introduces one 
new degree of freedom into the surface. We can manip- 
ulate this degree of freedom using any chain of X opera- 
tors connecting this face or "smooth defect" to a smooth 
boundary and any chain of Z operators encircling the 
smooth defect as shown in Fig. [6] We choose to call any 
such X chain Xl, and any Z ring Z^. This implies that, 
by definition, our logical qubit is initialized to |0l) as the 
surface is initially in the simultaneous +1 cigenstate of 
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FIG. 5: (Color online), a.) Locations in space and time, 
indicated by red dots, where and when the reported syndrome 
is different from that in the previous time step. Note that this 
is not a three-dimensional physical structure, just a three- 
dimensional classical data structure, b.) Optimal matching 
highly likely to lead to a significant reduction of the number 
of errors if bit-flips are applied to the spacelike edges. 

z z 



! prn 

6 o 6 o 




uJu 

i — o — 1 — o 



X 



FIG. 6: (Color online). Surface code with one additional 
degree of freedom introduced by not enforcing the stabilizer 
associated with one face (shaded). This face, or a region of 
such faces, is called a smooth defect. The degree of freedom 
can be phase-flipped by any ring of Z operators encircling the 
defect and bit-flipped by any chain of X operators connecting 
the defect to a smooth boundary. 



every Z stabilizer and therefore also in the +1 eigenstate 
of Z L . 

Larger smooth defects can be created using X mea- 
surements as shown in Fig. [7] Note that arbitrarily large 
defects still only introduce one degree of freedom. The 
given example shows the removal of four qubits and five 
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FIG. 7: (Color online). Surface code with one degree of free- 
dom introduced via the measurement Mx of four qubits in 
the X basis and removal of five stabilizers. Note that four 
new three term X stabilizers are created with not necessarily 
positive sign (indicated by green triangles). 



stabilizers — four Z stabilizers and one X stabilizer. Af- 
ter the X measurements, a number of new three term X 
stabilizers are created with not necessarily positive sign. 
As in the case of surface initialization, we will treat any 
negative eigenvalues as syndrome changes which will then 
be matched and corrected with chains of Z operators. 
The qubits inside the defect, which have been projected 
into a product state, play no further role in the compu- 
tation unless the defect moves. 

In practice, it is inconvenient to use a logical qubit with 
a logical operator that connects to a potentially distant 
boundary. This situation can be avoided by using a pair 
of defects to represent a single logical qubit as shown in 
Fig. [8l A chain of X operators connecting the two defects 
is then used as the Xl operator. The Zl operator is any 
ring of Z operators around either defect — these two 
classes of Zl operators are equivalent as they have the 
same commutation relations. 

Effectively, the above means we are choosing to repre- 
sent an arbitrary logical state by ct|0i) |0i) lx,) | lj,) as 
defined in the opening paragraph of this section. For the 
remainder of the review we shall redefine |0l) and 
such that an arbitrary logical state of a double defect 
logical qubit can be expressed as simply a|0i) + 
Note that double smooth defect logical qubits are also 
initialized to |0l) by default. 

Double smooth defect logical qubits can also be initial- 
ized in the state by first preparing a region of |+) 
as shown in Fig. [9] Such a region is automatically in the 
+1 eigenstate of Xl operators and X stabilizers not in- 
tersecting the boundary. X stabilizers on the boundary 
will have random sign. Smooth defects are then created 
by measuring all Z stabilizers outside the desired defect 
locations. The signs of the Z stabilizers will be random 
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FIG. 8: (Color online). Smooth qubit comprised of two 
smooth defects. Zl corresponds to any ring of Z operators 
around either defect. Xl corresponds to any chain of X op- 
erators connecting the two defects. 



FIG. 10: (Color online). Initializing a rough qubit in the 
state via Z basis measurements Mz and ignoring stabilizers 
(shaded). Xl is any ring of X operators around either defect. 
Zl is any chain of Z operators linking the two defects. 
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FIG. 9: (Color online). Initializing a smooth qubit in the \+l) 
state. After preparing a region of qubits each in the |+) state, 
every X stabilizer on the boundary of the region and every 
Z stabilizer outside the dashed regions is measured. Negative 
eigenvalues are treated as errors and corrected. 



and we will again treat negative stabilizers as syndrome 
changes, match them and correct them with chains of X 
operators. We will henceforth refer to a double smooth 
defect logical qubit as simply a smooth qubit. 

Rough qubits are also possible to create via Z measure- 
ments as shown in Fig. [TOl In this case the Zl operator 
is any chain of Z operators linking the two defects, and 
Xl any ring of X operators around either defect. Rough 
qubits are initialized to the +1 eigenstate of Xl, \+l), 
by default, although |0l) can be prepared starting with 
a region of qubits in the |0) state. 

Logical measurement is similar to initialization. To 
measure a smooth qubit in the Zl basis, a region of 
qubits encircling either or both defects is measured in the 
Z basis. In the absence of errors every path encircling ei- 
ther defect will have the same parity of Z measurements. 
If errors are present, they can be detected and corrected 
using the standard error correction procedure as directly 
measuring qubits in the Z basis is also an acceptable 
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FIG. 11: Example of measurement of a smooth qubit in the 
Zl basis in the absence of errors. Note that the measurements 
around every face have even parity whereas the parity of any 
path of measurements encircling either defect is odd. The 
figure thus corresponds to the measurement result |li). 



way to gain information about the eigenvalues of the Z 
stabilizers — even parity of Z measurements around a 
face corresponding to a positive eigenvalue and odd par- 
ity corresponding to a negative eigenvalue. Note that, as 
shown in Fig. 1111 it is possible for every face to have even 
parity, meaning no errors, and every path around either 
defect to have odd parity, meaning a readout result of 
|U>. 

A smooth qubit can be measured in the Xl basis by 
measuring a region including both defects in the X basis. 
In this instance the parity of all chains of X measure- 
ments connecting the two defects will be the same in the 
absence of errors. Similarly, rough qubits can be easily 
measured in either logical basis. 



V. LOGICAL CNOT 

So far, we have discussed two types of logical qubits, 
smooth and rough, schemes to initialize and measure 
them in the Zl and Xl bases, and Zl and Xl oper- 
ations. The only two logical qubit gate in this scheme 



6 



is the logical CNOT gate. To understand how logical 
CNOT works, we first need to understand in detail the 
effect of moving a smooth defect. 

Consider Fig. [T2k . This shows a smooth defect and 
two stabilizers — a single face Z stabilizer and a Zl sta- 
bilizer. If we now measure the center qubit in the X 
basis as shown in Fig. 112b. we will be left with the center 
qubit in the ±X eigenstate, and a stabilizer equal to the 
product of the face and the path. We have effectively de- 
formed the shape of the Zl stabilizer without changing 
its sign. The movement of the defect can be completed 
by measuring the Z stabilizer indicated in Fig. [P2b . and 
possibly correcting the sign of the result by applying an 
X operator to the center qubit. By repeating this pro- 
cess we can see that moving a smooth defect deforms the 
shape of Zl stabilizers passing nearby. 

Consider Fig. [T3k . This shows a smooth defect and 
three Xl stabilizers. Measuring the center qubit in the 

X basis as before, we see in Fig. H3b that this has po- 
tential side-effects, with a negative eigenvalue indicating 
the creation of three term negative X stabilizers and Xl 
stabilizers of changed sign. As shown in Fig. [T3h . the 
measured qubit, or qubits in the case of a larger defect, 
are individually phase-flipped to ensure they are all in 
the +1 eigenstate. Pairs of three term negative X stabi- 
lizers are corrected with chains of phase flips along the 
boundary of the defect. Note that this also corrects any 

XI stabilizers of changed sign. Fig. [T3H shows the effect 
of completing the movement of the defect by measuring 
the appropriate Z stabilizer. With the signs of the Xl 
stabilizers appropriately corrected, all Xl stabilizers at- 
tached to the defect remain attached to the defect with 
unchanged sign. By repeating this process we can see 
that moving a smooth defect drags around Xl stabiliz- 
ers attached to it. 

At first glance, the procedure described in the previ- 
ous paragraph does not appear to be fault-tolerant as it 
seems to rely on perfect measurement and correction of 
single qubits. Indeed, the procedure is not fault-tolerant 
unless a larger defect is used as shown in Fig. [Mb . After 
measuring a region of qubits in the X basis and using 
the individual measurements to give the sign of the X 
stabilizers across the entire measured region by taking 
their local parity, we use phase-flips to reset them to the 
+1 eigenstate as best as we are able. We do not assume 
that we achieve this perfectly. Resetting helps simplify 
the later incorporation of these qubits into the surface 
code. 

The three term negative X stabilizers on the boundary 
left over after resetting are again treated as syndrome 
changes and corrected using the procedure outlined in 
Section IIIII and described in more detail in Section IVIIII 
with the exception that if the procedure suggests con- 
necting to syndrome changes on the boundary, either new 
or old, the direction of the chain of operators correcting 
these changes is chosen such that the minimum number 
of sites on the old (reliable) boundary are changed. By 
doing this, every round error correction makes it expo- 
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FIG. 12: (Color online), a.) Smooth defect and surface in the 
+1 eigenstate of Zl- b.) After measuring the center qubit in 
the X basis, the shape of the Zl operator is deformed, c.) 
Measuring and possibly correcting the indicated Z stabilizer 
using a bit-flip on the center qubit completes the movement 
of the defect. 
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FIG. 13: (Color online), a.) Smooth defect and surface in the 
+1 eigenstate of Xl- b.) After measuring the center qubit in 
the X basis, it is possible that three term X stabilizers and Xl 
stabilizers with negative sign are created (potential locations 
indicated in green), c.) All signs can be corrected by applying 
the appropriate single qubit Z operators and chains of Z op- 
erators, d.) Measuring and possibly correcting the indicated 
Z stabilizer using a bit-flip on the center qubit completes the 
movement of the defect. 



nentially less likely that the three term negative X stabi- 
lizers created during the measurement step in a fixed unit 
length of boundary still remain. This implies that the 
number of rounds of error correction required to achieve 
a fixed probability of no errors remaining from the mea- 
surement step only grows logarithmically with the length 
of boundary. Note, however, that during the correction 
procedure new errors can occur. The primary desirable 
feature of these new errors is that they are unlikely to 
form very long chains. 

Fig. Q3t> shows a potential challenge when it is time 
to shrink the size of the defect and complete its move- 
ment. It is possible for a pair of errors to remain with 
one error on the boundary of the region of defect about 
to be healed, and the other error on the boundary of the 
region of defect to remain. As shown, without correction, 
this would result in Xl stabilizers with sign dependent on 
where they attach to the moved defect — a situation that 
is not allowed. However, by measuring all appropriate X 
stabilizers outside and on the boundary of the final posi- 
tion of the defect (indicated by a dashed red line) before 
measuring all of the Z stabilizers outside the final po- 
sition, the presence of these two errors is preserved and 
subsequent correction by a chain of Z operators ensures 
that the sign of all deformed Xl stabilizers is the same. 

To summarize the smooth defect movement procedure, 
a region of qubits is measured in the X basis, corrected 
as best as possible so that each measured qubit is in the 
+1 eigenstate, which takes constant time, several rounds 
of error correction are then applied until it is sufficiently 
likely that only new errors occurring after the initial mea- 
surements are now present on the boundary of the new 
region, which takes a time that grows logarithmically in 
the length of the boundary then measurement of all X 
stabilizers outside and on the boundary of the desired fi- 
nal position of the defect, measurement of all Z stabiliz- 
ers outside the final defect position, and finally error cor- 
rection until it is sufficiently likely that only errors occur- 
ring after the Z stabilizers were measured remain, which 
takes a time that grows logarithmically in the area being 
corrected. This movement procedure deforms nearby Zl 
stabilizers and drags around Xl stabilizers attached to 
the defect and takes a total time that grows only loga- 
rithmically in the distance the defect is moved. 

Now that we have a thorough understanding of the 
effect of moving a smooth defect, wc can return to the 
question of how to build a logical CNOT. Any gate can 
be completely specified by stating its action on computa- 
tional basis states, and can be specified up to global phase 
by stating its action on a basis of stabilizers. Specifically, 
if we have a system of two qubits and denote the CNOT 
between them with the first qubit as the control as A12, 
by simple matrix multiplication we can show that the 
following relationships hold 



A 12 (I®X)A\ 2 
A 12 {X®I)A\ 2 = 

K 12 {I®Z)K\ 2 : 



--I®X 

x®x 

--Z®Z 
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(7) 
(8) 
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FIG. 14: (Color online), a.) Movement of a large smooth defect via many measurements in the X basis. Many pairs of three 
term X stabilizers with negative sign are likely to be created (indicated in green), b.) After several rounds of error correction, it 
becomes exponentially unlikely that three term X stabilizers with negative sign remain that were generated in the measurement 
round. New chains of errors on the boundary can occur, but these will be corrected during normal error correction after the 
size of the defect is reduced to complete the movement. 



K 12 (Z®I)K\ 2 = Z®I (9) 

These relationships can be combined to determine the 
action of CNOT on an arbitrary two-qubit stabilizer. To 
show that we have a logical CNOT, it is sufficient to show 
that we can transform logical stabilizers in the above 
manner. Figs. fTSHTBl show that the logical versions of 
Eq. [7] and Eq. [S] hold if we use a smooth qubit as the 
control and a rough qubit as the target and braid one of 



the smooth defects around one of the rough defects. It 
is not important in which direction the braiding is done, 
only that the defect return to its initial position. It is 
also not important which smooth defect is moved nor 
which rough defect it is braided around. It is similarly 
straightforward to show that Eq. [5] and Eq. [5] hold. 

We do not yet have what we truly need — a CNOT be- 
tween logical qubits of the same type. Define Mx, Mz to 
take the value when the +1 eigenstate is measured and 
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FIG. 15: (Color online), a.) Surface containing a smooth qubit in the +1 eigenstate of Xl and a rough qubit. The lower 
smooth defect has been braided around the upper rough defect using X measurements. Note that is not possible to complete 
the braiding in one step as a ring of X measurements corresponds to measurement of the rough qubit in the Xl basis, b.) Via 
correction of many Z stabilizers, the Xl operator is dragged around the upper rough defect, c.) Additional X measurements 
extend the defect back to its original position, d.) Further correction of Z stabilizers returns the defects to their original 
positions but the surface is now in the +1 eigenstate of both the smooth and rough Xl operator. 



1 when the -1 eigenstate is measured. Consider Fig. [T7k . 
This is built entirely out of logical circuit elements de- 
scribed above and is equivalent to Z Mx on the target 
qubit followed by CNOT followed by X Mz on the target 
qubit. This is in turn equivalent to CNOT followed by 
(Z <S> Z) Mx followed by X Mz on the target qubit. We 



will adopt the policy of applying corrective logical op- 
erations based on the measurement results immediately 
after such a CNOT to simplify the discussion of more 
complicated circuits. Fig. [T7a can also be represented as 
a braiding of defects of different types in two dimensions 
of space and one dimension of time as shown in Fig. 117b 
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FIG. 16: (Color online), a.) Surface containing a smooth defect and a rough defect in the +1 eigenstate of Zl- The lower 
smooth defect has been braided around the upper rough defect using X measurements, deforming the shape of the rough Zl 
operator, b.) By first correcting many Z stabilizers and then performing further X measurements, the smooth defect can 
be extended back to its original position, c.) A final round of Z stabilizer correction returns the defects to their original 
configuration but with the state of the surface changed, d.) The Zl operator shown in part c is equivalent to the tensor 
product of smooth and rough Zl ■ 
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FIG. 17: (Color online). Smooth qubits are represented by 
black lines, rough qubits by red lines, a.) Circuit equiva- 
lent to Z Mx on the target qubit followed by CNOT between 
the control and target qubit followed by X Mz on the target 
qubit. b.) Schematic representing the initialization, braiding 
and measurement of defects in a surface code to implement 
Fig. 117b . Time runs from left to right, and the surface code 
should be imagined oriented vertically and into and out of the 
page, c.) Simplified schematic equivalent to Fig. 117b . 



and simplified in Fig. [17b. Note that since logical CNOT 
is built out of defect movement, which takes a time that 
grows logarithmically in the distance moved, and defect 
measurement, which takes constant time, the total time 
required to execute logical CNOT grows only logarithmi- 
cally in the separation of the logical qubits. 



VI. STATE INJECTION AND NON-CLIFFORD 
GATES 



The set of gates discussed so far is not universal. To 
complete the universal set, we will firstly describe how it 
is possible to non-fault-tolerantly prepare arbitrary logi- 
cal states, and then discuss state distillation 0, 0] and 
non-Clifford gates based on these distilled states. 



FIG. 18: Surface code fragment and numbered qubits used to 
assist the visualization of the discussion of Section [VI Al 



A. State injection 

Consider Fig. [TH We will focus on the numbered 
qubits and the four stabilizers X^XtX^Xq, 
Z2Z4Z5Z7, ZsZ^ZqZs centered on qubit 5. The discus- 
sion of this section applies to a surface of arbitrary size 
— we shall sec that none of the necessary manipulations 
affect stabilizers further away. We shall only explicitly 
work through the creation of an arbitrary rough qubit — 
the procedure for creating an arbitrary smooth qubit can 
be obtained by exchanging the roles of X and Z . 

To create an arbitrary rough qubit, begin by measuring 
qubit 5 in the X basis. This gives a state stabilized by 
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If the —1 eigenstate is obtained, apply either Z-iZ^Z^Zi 
or Z^Z^ZqZq to create the +1 eigenstate 
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Next, Hadamard transform (for pedagogical clarity) and 
then unitarily rotate qubit 5 to the desired state 
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If the —1 eigenstate of Z 2 Z i Z 5 Z- 7 and Z 3 Z 5 Z 6 Z 8 is 
obtained, apply X 5 and then either X1X2X3X5 or 
X^X-fX^Xg to give the desired logical state 
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After creating an arbitrary logical qubit using the pro- 
cedure above, the two halves of the logical qubit would be 
both moved apart and made larger as quickly as possible 
to make the logical qubit fault-tolerant. 



B. State distillation 

For our purposes, we are interested in the injection 
of two particular states \Y) = |0) + i\l) and \A) = 
|0) + e 171 "/ 4 1 1) . These states have very special properties. 
Consider Fig. [T9b. This takes an arbitrary state \ip) and 
six ancilla qubits initialized to |0) and creates a logical 
qubit protected by the 7-qubit Steane code [13] • A pe- 
culiar property of this encoding circuit is that if it is run 
backwards with seven states approximately equal to \Y) 
as in Fig. [T9b. the output will be closer to \Y). This 
is state distillation [l5|, [Ty| . Repeated multiple times, 
arbitrarily high fidelity \Y) states can be obtained ex- 
ponentially quickly. Specifically, if the input states each 
have a probability p of X, Y or Z error, the output state 
will have probability 7p 3 of error [2j • 

Some technicalities exist surrounding the actual val- 
ues of the measurements indicated in Fig. 119b . If perfect 
\Y) states are input, Table. U summarizes the possible 
measurement patterns, their probabilities, and the out- 
put state. It can be seen that in some cases a corrective 
Z operator needs to be applied to the output. If less than 
perfect |Y) states are input, other measurement patterns 
have nonzero probability. If a measurement pattern not 
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FIG. 19: a.) Encoding circuit for the 7-qubit Steane code, 
b.) Distillation circuit for the \Y) = |0) + state. 



listed in Table [I] is obtained, the distilled state is dis- 
carded. Note that for high fidelity input states the prob- 
ability of obtaining an unacceptable measurement pat- 
tern is asymptotically zero, thus few distilled states need 
to be discarded. 

A very similar distillation circuit exists for the \A) 
state. Fig. [20k shows the encoding circuit for the 15- 
qubit Reed-Muller code [l8j]. Running this backwards 
yields the distillation circuit shown in Fig. I2"0b . As be- 
fore, given perfect input states, only certain measurement 
patterns are possible and all measurement patterns result 
in the desired \A) state, though this time up to an X, Y 
or Z operator. Convergence is similarly rapid with error 
probability p input states yielding an error probability 
35p 3 output state [Hj]. 

Note that both Fig.fTUb and Fig.l2T)b are made of oper- 
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FIG. 20: a.) Encoding circuit for the 15-qubit Reed-Muller code, b.) Distillation circuit for the \A) = |0) +e" r/4 |l) state. 



Pr(M) Mx Mx Mx Mz Mz Mz |jg) 



0.125 




















Z\Y) 


0.125 








1 


1 


1 


1 


Z\Y) 


0.125 





1 





1 





1 


\Y) 


0.125 





1 


1 





1 





\Y) 


0.125 


1 











1 


1 


\Y) 


0.125 


1 





1 


1 








\Y) 


0.125 


1 


1 





1 


1 





Z\Y) 


0.125 


1 


1 


1 








1 


Z\Y) 



TABLE I: Possible measurement patterns after running the 
distillation circuit of Fig. 119b with perfect \Y) states and no 
gate errors. 



ations that can be performed easily and efficiently using 
the surface code. In particular, the single control multiple 
target CNOTs can be implemented in the same amount 
of time as a single CNOT. The input \Y) and \A) states 
would be created factory style, with any errors detected 
early in the non-fault-tolcrant process of their creation 
resulting in a restart of the creation process. Logical an- 
cilla states that are likely to be sufficiently good would 
then be recursively fed into logical surface code versions 
of Fig. [19b and Fig. [20b until sufficiently high fidelity 
ancilla states are obtained. 



C. Non- Clifford gates 

Given states of the form (|0) + e l6 \l))/y/2, rotations 
R z {6) and Rx{6) can be performed using the circuits 
shown in Fig. [2Tb and Fig. [2Tb respectively. Note that 
both of these circuits are probabilistic, and actually per- 
form rotations XRz(—0) and ZRx{—8) if the measure- 
ment indicates a negative eigenstate. If we wish to apply 
Rz(ir/2) or Rx(tt/2) and discover we have actually ap- 
plied rotations XRz(—ir/2) or ZRx{— t/2), the correct 
gate can be achieved simply by a subsequent application 
of Z and X. If attempting R z (n/4) and we discover we 
have applied XRz(— 7r/4), an ancilla state |0)+i|l) needs 
to be ready for an attempt to apply Rz(tt/2)X . If we 
again measure a negative eigenstate, subsequent applica- 
tion of ZX gives the desired rotation. 



VII. LOGICAL HADAMARD 

The Hadamard gate is called for in many quantum 
algorithms. In principle we could simply use the relation 

H = R z (-k/4)R x (-k/4)R z {-k/4) (15) 

and the constructions of Section IVTI There is, however, 
a much more efficient way (l2l |. 
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FIG. 21: a.) Circuit performing the single qubit 
unitary X Mz Rz({— l) Mz 0) given an appropriate ancilla 
state. b.) Circuit performing the single qubit unitary 
Z Mx Rx{{ — 1) x ff) given an appropriate ancilla state. 



Consider Fig. [22] This shows a smooth qubit cut out 
of a larger lattice using Z measurements. Note that the 
three term Z stabilizers thus created would need to be 
corrected as they would have random sign. Without cor- 
rection, the indicated Zl stabilizer would have random 
sign after the measurements. Note that the ring of Z 
measurements provides no information about the state 
of the smooth qubit — such a ring is equivalent to the 
logical identity operator. 

The logical Hadamard gate can now be performed 
transversely. Every face Z stabilizer becomes a vertex 
X stabilizer. The rough boundary becomes a smooth 
boundary. The smooth qubit becomes a rough qubit. 
Stabilizers Zl and Xl are interchanged. This last point 
is precisely the action of logical Hadamard. 

The interchanging of faces and vertices does create a 
slight problem — faces and vertices are no longer where 
they should be. Before connecting the logical qubit to the 
rest of the lattice, it would need to be moved diagonally 
in any direction a half lattice spacing. This could be 
achieved via physical swap gates. After realignment, the 
complete surface of stabilizers would be measured and 
corrected once more. 

We are still not quite done — our logical qubit has 
been converted from smooth to rough. By preparing 
a smooth ancilla qubit in the state and perform- 

ing a simple smooth-rough CNOT followed by measure- 
ment of the rough qubit in the Zl basis and application 
of Xl if the —1 eigenstate is obtained, we can convert 
the rough qubit back into a smooth qubit, completing 
the process. While not completely trivial, this complete 
process is vastly simpler than the necessary ancilla state 
preparation and distillation associated with Eq. (|15[) . 



VIII. THRESHOLD ERROR RATE 

In our simulations we look at a planar square lattice 
with two smooth and two rough boundaries. This type 
of a lattice lets us encode one logical qubit. Our general 
calculation strategy involves preparing the system in the 
simultaneous +1 eigenvalue of all Z and X stabilizers 
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FIG. 22: (Color online). A smooth qubit isolated from a 
larger piece of surface code using a ring of Z measurements 
so that logical Hadamard can be applied by local, transver- 
sal Hadamard gates. Red triangles represent three term Z 
stabilizers of negative sign. See text for details. 



and observing how long it takes for the encoded logical 
state to change as a result of randomly generating errors. 

We have not yet discussed how X and Z stabilizers are 
actually measured. Fig. [23] shows that a fifth syndrome 
qubit is required to detect whether the state of the sur- 
face \ip) is in the ±1 eigenstate of a Z or X stabilizer 
[l9| . If the surface is in neither eigenstate, the circuit 
projects the surface into a state \tp') that is one of the 
±1 eigenstates. It takes six steps to perform such a mea- 
surement. The syndromes are initialized, a CNOT op- 
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FIG. 23: Circuit showing how an additional syndrome qubit 
(top line of each figure) is used to measure a.) Z stabilizers, 
b.) X stabilizers. 



eration is applied between each syndrome and the qubit 
to the north, west, east and south, and finally every syn- 
drome qubit is itself measured, as schematically shown 
in Fig. [24] This order of the CNOT gates has been cho- 
sen to ensure that adjacent syndrome circuits sharing a 
pair of data qubits are strictly ordered — one syndrome 
circuit touches both data qubits before the other syn- 
drome circuit. CNOT gate orders without this property 
resulted in entangled syndrome qubits that provide no 
useful information for error correction. The placement of 
a syndrome qubit on each vertex and in the center of each 
face, plus the CNOTs required during the circuit imply 
that we need a 2-D nearest neighbor coupled lattice of 
qubits. 

The threshold error rate is derived from four error rates 
in our simulations — initialization error p^, readout error 
p r , memory error p m and the error associated with a two- 
qubit gate p g . Note that we combine any single-qubit 
gates with neighboring two-qubit gates and thus do not 
have a separate single-qubit error rate. All four of these 
error rates are set to the same value p and all operations 
are assumed to take the same amount of time to permit 
our threshold error rate to be compared with others in 
the literature [lEIUl. 

By initialization, we mean initialization to the state 
|0). An initialization error is therefore accidental prepa- 
ration of state |1) with probability p. By readout, we 
mean readout in the Z basis. A readout error is a classi- 
cal error — the qubit is projected into the ±1 eigenstate 
of Z, but with probability p the eigenstate reported by 
the measurement device is incorrect. A memory error is 
the application of X, Y or Z, each with probability p/3, 
to an idle qubit. A two-qubit gate error is the application 
of one of the 15 nontrivial tensor products of /, X 1 Y and 
Z, each with probability p/15, after perfect application 
of the two-qubit gate. 

As was briefly outlined in Section ITTTl after each syn- 




FIG. 24: (Color online). Syndrome measurement typically 
involves six gates: syndrome initialization, CNOTs with the 
four surrounding data qubits (fewer on boundaries) and finally 
syndrome readout. 



drome is read, its value is checked against a result from 
the previous iteration, and if the values differ, the syn- 
drome change location (in time and space) is recorded. 
Next, a matching of all the syndrome changes collected 
up to this point (an example is shown in Fig [5^) is used 
to guess where errors occurred. Since shorter error chains 
are more likely than longer ones, we use a minimum 
weight matching algorithm to do this [3]. Before the 
matching algorithm can find a minimum weight solution, 
we convert all the syndrome change results into a graph, 
with locations of the syndrome changes representing the 
nodes, and edges between these nodes having a weight 
which depends on the distance between them. The edge 
weight is measured in faces along the special dimensions 
and syndrome extraction cycles along the time dimen- 
sion. 

We once again stress that some error chains may begin 
at the boundary and end somewhere inside the lattice 
(see Fig [3]). In such cases, we can only observe the syn- 
drome change on the interior of the lattice. To account 
for this (meaning enable the matching algorithm to guess 
that the error chain started on a boundary), for every in- 
terior node, we always create a closest boundary node. 
The edges between different boundary nodes are set to 
be of weight zero. One way to prepare our graph would 
be to include an edge between every pair of nodes (since 
in principle we don't know where actual errors occurred), 
but in practice this is not necessary. Only edges that con- 
nect nodes which are not further from each other than the 
sum of the weights between each node and their closest 
boundary nodes are included since nodes that are further 
apart will always be matched with their respective closest 
boundaries in preference to each other. 

We further optimize graph creation by noting that 
matches which are temporally far behind the current 
time step are unlikely to be modified by recent syndrome 
changes and therefore can be "remembered" from pre- 
vious iterations. These techniques let us minimize the 
size of the graph that is passed to the matching algo- 
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FIG. 25: (Color online). Logical Z (X) error detection in- 
volves checking if the parity of Z (X) operators along any 
of the vertical (horizontal) lines of qubits is odd. Above, we 
show this in an example of a logical Z error. 



FIG. 26: (Color online). A plot of average time until failure 
versus the physical error rate p. A threshold is observed at 
p ~ 6.0 x 10 -3 where the curves of different lattice sizes cross. 
The case where no error correction is used (single qubit) is 
represented by a dashed line. 



rithm which, despite scaling polynomially in the number 
of edges and nodes, can often still take substantial com- 
puting time. An example of a successful minimum weight 
perfect match is shown in FigOj. 

As outlined at the beginning of this section, in order 
to know if the simulation should continue or not, we need 
to determine whether the lattice suffered a logical error 
(and hence the encoded state has changed). A logical 
error corresponds to a chain of errors that starts on one 
boundary and ends on the opposite one. In order to 
detect if a logical error has occurred, we repeat the read- 
out cycle with all the error sources set to zero (i.e. set 
Pi = Pr = Pm = Pg = 0, in other words have a "perfect 
readout"). This allows us to be certain that any logical 
Z (X) error can be recognized by solely checking if the 
parity of Z (X) operators crossing a vertical (horizontal) 
line of qubits is odd. A simple example of this is shown 
in Fig [25j If no logical error is detected, we revert the 
simulation state to what it was just before the "perfect 
readout" cycle was executed and continue on. 

During every run, we note how many syndrome ex- 
traction cycles it took for a logical error to be observed. 
The simulation is repeated many times for different lat- 
tice sizes and values of the physical error rate p. All this 
data is then used to calculate the average number of steps 
until a logical error occurs for a given lattice size and p. 
The graph of Fig |2"6"1 shows the obtained results. In it we 
see a log-log plot of the average time until failure versus 
the physical error rate p for lattice sizes ranging from 
4 to 20 faces across. We observe a crossing at approxi- 
mately p«6.0xl0~ 3 , which is our numerical threshold. 
If the physical error rate is below this threshold value, 
the average number of readout cycles until failure can 
be increased arbitrarily by increasing the distance of the 
code (lattice size). 
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FIG. 27: (Color online). Potential X error chains (dashed 
lines) around a rough defect and consequent minimum sepa- 
ration from long straight smooth boundaries and two types of 
corners such that X error chains beginning and ending on the 
smooth boundary are not more likely than an X error ring 
around the rough defect. 



IX. DISTRIBUTED COMPUTING 

In this section, we show that distributed quantum com- 
puting can be performed in a natural manner. For our 
purposes, a distributed quantum computer will consist of 
a number of separate rectangular lattices of qubits each 
capable of holding at least two logical qubits. Comput- 
ing shall proceed by first moving logical qubits that need 
to interact onto a common plate before attempting the 
logical interaction. The movement of logical qubits from 
one plate to another is the only additional capability we 
need to discuss. 

Consider Fig. [28k . This shows a plate containing a 
rough qubit and an empty plate. Note that rough de- 
fects do not need to be kept very well separated from 
smooth boundaries as no error chain can link a rough de- 
fect with a smooth boundary. Fig.[57]shows the minimum 
permissible separation from long straight boundaries and 
corners. Rough defects do, however, still need to be kept 
well separated from each other. 

To move the rough qubit from one plate to the other, it 
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must be possible to perform remote gates between either 
two complete edges, or a smaller section of two edges 
if the plates are large relative to the size of a logical 
qubit. Generally speaking, implementing remote gates 
would be expected to involve entanglement distribution 
and purification [2l[ . We will not discuss the details here 
besides mentioning that this leads to significant qubit 
and gate overhead implying remote gates should be kept 
to a minimum. 

Consider Fig. I2"5b . This shows the pairs of qubits, 
including syndrome qubits, that need to be remotely in- 
teracted to enable a single round of the error correction 
to proceed seamlessly across the two plates. Note that 
one column of qubits on the empty plate has been omit- 
ted as though it is idle, but note that this figure does 
not include the qubits required for entanglement purifi- 
cation and it is unlikely there would be idle qubits on the 
boundary in practice. In general, the joined plates will 
be in random eigenstate of both the X and Z stabiliz- 
ers straddling both plates. We shall treat these random 
values as errors and correct them. 

After correction of the join, the rough qubit can be 
moved over to the other plate via Z measurements as 
shown in Fig. [28b . First the border Z stabilizers of this 
extended defect would need to be corrected as discussed 
in Section [Vj then, when shrinking the size of the de- 
fects to move the logical qubit, the unneeded regions of 
X stabilizers measured and corrected once more. Both of 
these correction procedures take a number of time steps 
that only grows logarithmically with the size of the com- 
putation and the length or area being corrected. After 
the necessary correction has been completed, error cor- 
rection can continue on each plate individually without 



any further long-range interactions. 

The most common reason to move a logical qubit 
from one plate to another would be to perform a remote 
CNOT. This would be achieved by creating a rough qubit 
on the control plate, braiding it around the control qubit, 
sending the rough qubit to the target plate and complet- 
ing the necessary braiding and measurement operations 
entirely on the target plate. 



X. CONCLUSION AND FURTHER READING 

We have presented a simplified yet comprehensive re- 
view of the 2-D version of the quantum computation 
scheme originally presented in [H, 0]. We started with 
a description of the surface code, as well as the stabilizer 
formalism which is used throughout this paper. We dis- 
cussed in detail logical state initialization, logical CNOT 
and non-Clifford group gates, which make use of state 
distillation. We calculated a numerical threshold for the 
surface code and obtained a value of p w 6.0 x 1CP 3 which 
is commensurate with other calculations in the literature 

atom. 
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FIG. 28: a.) A rough qubit ready to be sent to a separate piece of surface code, b.) Remote gates are used to join the two 
surfaces together, c.) A sequence of measurements is used to move the rough qubit. After the necessary correction associated 
with completing the movement, the long-range gates can be discontinued to separate the two pieces of surface code once more. 



